Switching dynamics of surface stabilized ferroelectric liquid 
crystal cells: effects of anchoring energy asymmetry 

Alexei D. Kiselev 

Institute of Physics of National Academy of Sciences of Ukraine, 
prospekt Nauki ^6, 03028 Kyi'v, Ukrain^ 

Vladimir G. Chigrinov 

Hong Kong University of Science and Technology, 
Clear Water Bay, Kowloon, Hong Kong 

Eugene P. Pozhidaev 

P.N. Lebedev Physics Institute of Russian Academy of Sciences, 
Leninsky prospect 53, 117924 Moscow, Russia 
(Dated: February 6, 2008) 

We study both theoretically and experimentally switching dynamics in surface 
stabilized ferroelectric liquid crystal (SSFLC) cells with asymmetric boundary con- 
ditions. In these cells the bounding surfaces are treated differently to produce asym- 
metry in their anchoring properties. Our electro-optic measurements of the switching 
voltage thresholds, V+ and — V-, that are determined by the peaks of the reversal po- 
larization current reveal the frequency dependent shift of the hysteresis loop, V+ — V-. 
We examine the predictions of the uniform dynamical model with the anchoring en- 
ergy taken into account. It is found that the asymmetry effects are dominated by 
the polar contribution to the anchoring energy. Frequency dependence of the voltage 
thresholds is studied by analyzing the properties of time-periodic solutions to the 
dynamical equation (cycles). For this purpose, we apply the method that uses the 
parameterized half-period mappings for the approximate model and relate the cycles 
to the fixed points of the composition of two half-period mappings. The cycles are 
found to be unstable and can only be formed when the driving frequency is lower 
than its critical value. The polar anchoring parameter is estimated by making a 
comparison between the results of modelling and the experimental data for the shift 
vs frequency curve. For a double-well potential considered as a deformation of the 
Rapini-Papoular potential, the branch of stable cycles emerges in the low frequency 
region separated by the gap from the high frequency interval for unstable cycles. 
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I. INTRODUCTION 



As it was shown by Meyer at al. in 1975 the origin of ferroelectric ordering in smectic- 
C* (Sm-C*) liquid crystals with the spontaneous polarization normal to the tilt plane is 
closely related to the reduction of symmetry caused by chirality of the tilted Sm-C* phase. 
This phase is thus naturally ferroelectric and the Sm-C* liquid crystals are also known as 
the ferroelectric liquid crystals (FLC) (a detailed description of FLCs can be found e.g. in 
a collection of reviews [2] and in a more recent monograph [3]). 

Over the past three decades FLCs have been attracted considerable attention not only as 
a vital issue in the condensed matter physics but also as promising materials for applications 
in electro-optic switching devices. The first such device was due to Clark and Lagerwall 0| ■ 
They studied electro-optic response of FLC cells confined between two parallel plates subject 
to homogeneous boundary conditions and made thin enough to suppress the bulk chiral 
Sm-C* helix. It was found that such cells — the so-called surface- stabilized ferroelectric 
liquid- crystal (SSFLC) cells — exhibit high-speed, bistable electro-optical switching between 
orientational states stabilized by surface interactions. 

The response of chiral Sm-C* liquid crystals to an applied electric field is characterized 
by fast switching times due to linear coupling between the field and the spontaneous polar- 
ization. There is also a threshold voltage necessary for switching to occur and the process 
of bistable switching is typically accompanied by a hysteresis. The switching times and the 
threshold voltages may considerably vary depending on the wave form, the amplitude and 
the frequency of applied (chiving) voltages. 

In early studies 0,0,0], it was found that the competition between elastic, electrostatic 
and surface energies may result in different regimes of switching and field induced transitions 
in FLC cells. Certain regimes such as the high voltage regime can be described using 
the theoretical approach based on the assumption that director and polarization fields are 
spatially homogeneous. This approach provides uniform switching models where the effects 
of electrostatic and surface interactions are incorporated into an effective potential governing 
the dynamics of switching. Such models can be readily applied to interpret experimental 
data. 

In particular, the uniform model of switching supplemented with an elastic-like term 
was examined in as a useful tool to describe experimental behavior of the polarization 
reversal current in SSFLC cells. Similarly, the uniform theory was applied to determine 
the rotational viscosity and the anchoring energy strength from the experimental results 
on the response time measured as a function of pulse voltage in SSFLC cell |9[. The field- 
reversal method suggested in [10| to measure the spontaneous polarization, the switching 
time, the rotational viscosity and the dc conductivity also relies on the uniform model. 
(Recent discussions, applications and generalizations of uniform models can be found e.g. 

In this paper we are aimed to study the effects of surface anchoring energy in switching 
dynamics of asymmetric SSFLC cells where nonidentical aligning films impose different 
boundary conditions at the substrates. The effect of our particular concern is the frequency 
dependent shift of the hysteresis loop observed in our experiments. 

In the theoretical part of the paper, we adapt a systematic approach and examine pre- 
dictions of the uniform theory where the effects of asymmetry are caused by the polar 
contribution to the anchoring energy potential. The voltage thresholds in relation to the 
driving frequency are studied by using the method suggested to explore the properties of 



periodic solutions to the dynamical equation. 

The layout of the paper is as follows. In Sec. HU we derive the effective potential and 
formulate the model. The simplified case with the anchoring energy neglected is discussed 
so as to clarify the assumptions taken to obtain frequency dependent threshold voltages. 

Analysis of the switching dynamics is performed in Sec. IIII1 We present the approach 
that uses the parameterized half-period mappings for the approximate model to study the 
time-periodic solutions of the dynamical equation. These solutions are called the cycles and 
related to the fixed points of the composition of the half-period mappings. The method 
is applied to the case of square wave voltages. Analytical relations for the conditions of 
complete switching, the switching times and the critical frequency bounding the region of 
cycles from above are obtained. For sine-wave and triangular voltages, we comment on the 
numerical procedure and show that, qualitatively, the results remain unchanged. 

Experimental details and the results of electro-optical measurements for the switching 
voltage thresholds are given in Sec. IIVI The experimental data and the theoretical results 
of Sec. IIHI are used to model the process of switching and to estimate the polar anchoring 
parameter. Finally, in Sec. 0we present our results and make some concluding remarks. In 
the Appendix we extend our analysis to the dynamical model with the double-well effective 
potential and study the branches of stable and unstable cycles. 

II. MODEL 

In this section we introduce necessary notations and describe the model that takes into 
account the effects due to aligning films. The equation of motion governing switching dynam- 
ics in asymmetric SSFLC cells are derived and the asymmetry induced polar contribution to 
the anchoring energy is found to play an important part in the problem Frequency depen- 
dence of the switching voltage thresholds is analyzed using a simplified model as an example 
demonstrating the key points underlying subsequent theoretical considerations. 

A. Free energy 

We consider the FLC cell with a planar arrangement of planar smectic layers usually 
referred to as the bookshelf geometry. This geometry is schematically presented in Fig. [H 
It is characterized by the director field 

n = z cos 8 + sin8 (x cos<ft + y sin0), k = y, (1) 

where 9 is the molecular cone angle, <fi is the azimuthal angle around the smectic cone, k is 
the outward (inward) normal to the upper (lower) substrate of FLC cell, y = d/2 (y = —d/2), 
and d is the cell thickness. 

The vector of the polarization, P, and the electric field inside the cell, E, are given by 

P = P s p, p = y cos0 - x sin0, E = Py, (2) 

where Ps is the spontaneous ferroelectric polarization. The external driving voltage 



V(t)=V m v(t) 



(3) 



FIG. 1: Schematic representation of bookshelf geometry illustrating planar arrangement of smectic 
layers. The polarization vector lies in the x — y plane and <f> is the azimuthal angle around the 
smectic cone. External electric field is directed along the y-axis normal to the boundary surfaces, 
y = —d/2 and y = d/2. 

is characterized by the amplitude V m , the frequency / = 1/T (T is the period) and the 
T-periodic function of time v(t) describing the waveform of the applied voltage. 
For sinusoidal (sine- wave) driving voltages, the waveform function v(t) is given by 



v S me(t) = sin(W), lu = 2n/T, 



(4) 



whereas T-periodic continuations of the functions ftriang(^) and v. 



squar 



(t): 




(5) 



(6) 




give the waveform functions for the triangular and square-wave forms, respectively. 

Similar to nematic liquid crystals (NLC), when a FLC layer is brought into contact with 
aligning substrates, the energy of the FLC molecules in the interfacial layer and thus the 
surface tension will be orient at ionally dependent. The anisotropic part of the surface ten- 
sion — the so-called anchoring energy — gives rise to the phenomenon known as anchoring, 
i.e., surface induced alignment of the FLC director along the vector of preferential orientation 
referred to as the easy axis. 

Phenomenologically, the expression for the anchoring energy can be written as a linear 
combination of the invariants constructed from the surface normal, k, the FLC director, n 
and the unit polarization vector p. For FLC cells, this gives the anchoring energy expressed 
in terms of the scalar products (n • k) and (p • k) as follows 



anch 



v=±l 



^(n.k) 2 + ^f)(p-k) 



(?) 



. y=vd/2 



where W + (W_) is the strength of nonpolar anchoring at the upper (lower) substrate; sim- 



ilarly, and W^f ' is the polar anchoring strength at the upper and lower substrate, 

respectively. 

The first term in square brackets on the right hand side of Eq. (jZJ) represents the Rapini- 
Papoular surface potential preserving equivalence between n and — n. This equivalence, 
however, can be broken due to effects of polar ordering in the interfacial layer Rj 1^ - 
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FIG. 3: Hysteresis loop. 

As opposed to the case of nematic liquid crystals, where the effects of surface induced 
polarity are mainly caused by the quadrupole-dipole interaction [l3,E^, one might expect 
that in FLCs the dominating factor is the electrostatic interaction of the spontaneous po- 
larization P with the surface charges or dipoles. So, the polar anchoring term on the right 
hand side of Eq. (JJJ) is taken to be proportional to (p ■ k) . This term introduces dependence 
on the polarity of the polarization resulting from the ferroelectric polar surface interaction. 

When the director field (JTJ is spatially homogeneous, the anchoring energy (J7J) takes the 
simplified form 

W 



W, 



anch 



sin — Wp cos < 



(8) 



where W = (W_ + W+) sin 2 9 and W P = - W^K From Eq. © it is seen that, for 



r(P) 



symmetric cells with W. 



(P) 



(P) 



the polar terms add to zero. In asymmetric FLC cells, 



the polar anchoring parameter, w p = Wp/W, is generally non- vanishing. 

In Figure the potential (JEJ) is plotted in relation to the azimuthal angle 0. It is 
illustrated that there are two local minima: = (the up state with p = (p ■ k) = 
cos0 = +1) and = tt (the down state with p = —1). When w p ^ 0, the minima are not 
energetically equivalent and one of the minima represents the metastable state. In addition, 
it is not difficult to obtain the relation 



\Wr. 



\W P /W\ < 1 



(9) 



as the condition for the configurations with = and = 7r to be metastable. If the 
metastability condition Q is broken the anchoring potential has only one minimum. Similar 
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FIG. 4: Normalized polarization and polarization reversal current evolving in time under sine- wave 
driving voltage at w p = 0.2, T w = 0.8 and r e = 50. Solid lines represent the results obtained by 
solving Eq. (|21j) numerically. Analytical results (dashed lines) are computed using the piecewise 
linear approximation (|40|) . 



results were previously reported in Refs. nana and the effects of asymmetry induced by 
the polar anchoring in switching dynamics will be of our primary concern. Throughout 
the paper we shall assume that, in the zero-field case, the up and down states are both 
metastable and the condition @ is satisfied. 

In general, the polar anchoring energy is known to be of considerable importance in the 
understanding of static and dynamical behavior of chiral smectics in confined geometries. 
For example, it may lead to the discontinuous Freedericksz transitions and is found to 
give rise to the surface electroclinic effect (2ll . 22l |. The magnitude of the polar anchoring 



strength was also estimated from the experimental data on the static dielectric susceptibility, 
the voltage coercivity and the relaxation tim e | 23l| . 

In a more recent theoretical investigation [24| into the effects of polar anchoring, which 
induces ferroelectric ordering close to the bounding surfaces of an antiferroelectric liquid 
crystal cell, it was found that such ordering may result in a coexistence of multiple zero- 
voltage ground states. 

As far as the SSFLC cells are concerned, the polar part of the anchoring energy was 
considered theoretically in the papers [13, HE, dealing with the so-called "thresholdless" 
or "V-shaped" switching mode of the optical response to an applied voltage. It was found 
that monostable structures exhibiting the V-shaped switching can be formed as a result of the 



polarization charge self-interaction accompanied by the effects of the insulating alignment 
layers. 

To see this, we briefly discuss a simple double-layer model. The starting point is the 
standard expression for the electrostatic part of the FLC free energy per unit area 



fd/2 

F E /A= / f E dy, f E = -e e yy E 2 /2-P y E, (10) 

J-d/2 



where A is the area of the substrates, Eq = 8.854 x 10 -12 Fm" 1 is the permittivity of free 
space and = e±5ij + (ey — e±)nirij is the FLC dielectric tensor. 

For FLC sandwiched between two insulating layers, we neglect the conductivity of FLC 
and use the boundary conditions for the normal component of the electric displacement field, 



D y , to obtain the relations 



e e yy E + P y = SqEiEi = e e 2 E 2 = D y , (11) 

where E\ {e<i) is the dielectric constant of the lower (upper) insulating layer and E% is the 
electric field inside the layers. The voltage applied across the cell is given by 



fd/2 

V = V(-d/2-d 1 )-V(d/2 + d 2 ) = E 1 d 1 + E 2 d 2 + Edy, (12) 

J-d/2 



where di is the thickness of the layers. 

We can now use Eq. (fTTj) to eliminate the electric fields, E\ and E 2 , from the expres- 
sion (fT2*|) and relate the electric displacement field to the external electric field, Eq, as follows 

V/d = E = V D y /E + (E) , V = [d x /e x + d 2 /E 2 )/d, (13) 

where (...) = d^ 1 f_ji 2 ■ ■ ■ dy. This result can be combined with the relation (fTTjl to yield 
the expression for the electric field inside the FLC layer 

h = FT — iT~i — i • ^ 14 ^ 

For a uniform director distribution, Eq. ([Tj| can be simplified giving the relation 

E= E - V P y /E 

1 + EyyT] 

recently obtained in Ref. j2?| for the case of identical substrates. Evidently, this relation 
implies that the electric field inside the FLC layer, E, deviates from the applied electric field, 
Eq, due to the presence of the aligning substrates. As it was emphasized in Refs. 0, 128^ . 
this effect plays an important part in describing the V-shaped switching mode. 

Alternatively, Eq. (fTKj) can be used to introduce the effective electrostatic free energy 

2/eff = -(1 + r,E yy )- l [EQE yy El + 2E P y - vP y /e ] (16) 

which is defined so as to meet the condition 



(17) 



expressing equivalence of the torques computed from the free energy (fTUj) and the effective 
energy (jl6j) . The last term in the square brackets on the right hand side of Eq. (JTfij) represents 
the energy of self-interacting polarization charges. 

Now we neglect the dielectric anisotropy of FLC, e yy ~ £flc = £± 5 and consider the 
energy 

F = fd = f e ffd + PVanch (18) 

defined as a sum of the effective and the anchoring energies. After substituting Eqs. (jHJ) 
and (fTfijl into the energy (JTHJ), it can be seen that there are two effects caused by the 
aligning layers: (a) reduction of the applied electric field: E — > [1 + t^flc] -1 -^ and 
(b) renormalization of the nonpolar anchoring strength induced by the polarization self- 
interaction: W — > W — 77 [1 + ^flc] l£ o -^i- The latter is the effect that may result 
in breaking the metastability condition (jHJ) thus leading to the formation of monostable 
structures exhibiting the thresholdless switching 113,12^. Since such a possibility as well as 
V-shaped switching is beyond the scope of this paper, we shall use the energy (fTHjl without 
changing notations for the renormalized values of Eq and W. 

The simple model of switching dynamics in FLC cells can be formulated in terms of the 
energy (JTHj) which enters the dynamical equation for the azimuthal angle 

90 df 

1 m=~W (19) 

where 7 is the rotational viscosity for reorientation on the smectic cone. The model is based 
on the assumption of spatially uniform director field that implies neglecting elasticity effects 
and the coupling between the fluid flow and the director. 

Such uniform theory suggests applying the viscosity-limited dynamics to the switching 
process in SSFLC cells. The simplest case where the anchoring energy is disregarded was 
originally considered in 0, H| (see e.g. for a review) . An attempt to fit experimental 
data using this simplified model was recently made in |l4(. 

Subsequently, various modifications of the model such as including an effective elastic-like 
term [H and the anchoring energy jj| 0, |26| have been employed to interpret experimental 
results. The approach to grey levels in FLC displays based on the uniform theory is presented 
in 0. 

In our case there are two characteristic time scales 

and the governing equation (fTTIJ) can be written in the following explicit form 

— = -(r(r) + cos + w p ) sin 0, (21) 
r = t/t w , r(r) = r e v(r), (22) 

where r e = tw/tE = PsVm/W is the driving voltage parameter. Note that rescaling of time 
renormalizes the period of the waveform function v(t): T —>■ T w = T/tw- 

Eq. (f2~Tj) presents the simplest model that can be used to study the effects of asymmetry 
induced by the polar anchoring energy. The characteristics of the hysteresis loop depicted 
in Fig. El will be of our primary concern. 



Referring to Fig. 01 location of the peaks of the polarization reversal current is deter- 
mined by two switching voltage thresholds, V+ = V m v + and V- = — V m V-, where v± is the 
dimensionless switching voltage parameter, characterizing the hysteresis loop for the nor- 
malized polarization (polarization parameter), p = P y /Ps — cos (p. Equivalently, width and 
shift of the loop can be conveniently described by the dimensionless (normalized) voltage 
parameters 

v c = v + + v-, v Bh = v + -V-, (23) 

where v c is the voltage coercitivity and t> s h is the voltage shift. 

If w p = 0, the dynamical equation (}2~T]) is invariant under the symmetry transformation: 
r — > — r (r — > r + T w /2) and <fi — > tc — (ft (p — > —p). Under these circumstances, the hysteresis 
loop is symmetric and the voltage shift vanishes, v s ^ = 0. This is no longer the case when 
w p ^ and the symmetry is broken. The curves representing temporal evolution of the 
normalized polarization, p, and its derivative with respect to time driven by the sine-wave 
voltage at w p = 0.2 are shown in Fig. |3J 

B. Switching voltages vs frequency: simplified model 

It is instructive to examine first predictions of the extremely simple model in which the 
terms describing the anchoring energy (jHJ) are disregarded and dynamics of the azimuthal 
angle <ft is governed by the simplified equation of motion 

l^ = -PsE v(t)sm<J ) . (24) 

Our task is to study the switching voltages as functions of the driving voltage frequency, 
/ = 1/T. We shall perform analysis of the model (J2*4"j) in the form suitable for subsequent 
generalization and clarify the assumptions underlying the results previously reported in 
Refs. 0,Ej. 

We begin with introducing an auxiliary angular variable u linked to the azimuthal angle 
through the relations 

2u = In ( — — ) , cos0 = tanh(w) (25) 

\ 1 — COS (p J 

and governed by the equation taken in the dimensionless form 

du , . . 

q^ = v(t), r = t/t E . (26) 

Solving Eq. gives the normal component of the polarization vector, P y , evolving in time 
as follows 

P y (r) = P s cos0(r) = P s tanh(«(r)), (27) 

u(t) = -u + w(t) = -u + / v(r')dr', (28) 

Jo 



where the initial condition cos (f)(0) = — tanh uq with uq > means that the cell is initially 
in the down state. 



Periodicity of the polarization ()27|) in time requires the condition w(T e ) = 0, where T e = 
T/tfi is the renormalized period, to be fulfilled. For bipolar switching with the waveform 
functions (jlj)- (0), we have 

w sine (r) = [1 - cos(^ e r)]/^ e , u e = 2vr/T e , T e = T/t E , (29) 

f 2(r/T e ) 2 , < r/T e < 1/4, 

^triang(r) = T e ^ -2(r/T e ) 2 + 2r/T e - 1/4, 1/4 < r/T e < 3/4, (30) 
[2(r/T e -l) 2 , 3/4 < r/T e <1, 

J r / T ^ 0<r/T e <l/2, 
^ squar(Tj " e \l-r/T e , 1/2 < r/T e <1, ( j 

so that the periodicity condition is satisfied. 

From Eq. (|2*Hj) the angular variable u monotonically increases from — uq to u + = —Uq + 
w(T e /2) over the first half-period. Complete switching occurs when cos0o — — tanhw ~ — 1 
and cos0 + = tanh-u + « 1. In this case, the polarization increases (decreases) passing zero 
at the instant of time r + (r_). So, the switching voltage parameters, v + and f_, can be 
defined as follows 

v± = ±v(t±), w(r±) = 0. (32) 

For sine- wave and triangular voltages, Eq. (|3*2*|l can be combined with Eqs. (|2*5|) -(|3*n |) to yield 
the expressions for the switching voltages: 

v± ne) = ^u u e (2 - u uj e ) } u uj e < 1, (33) 



(triang) _ \ a/8Mo/T 6 , < U /T e < 1/8, 

V± ~ \y/2 -8u /T e , 1/8 < u /T e < 1/4. 

It would appear natural that the polarization in the up and down states is of the same 
magnitude but differ in sign. So, the parameter u entering the expressions (|3*3*|) and 
can be fixed by imposing the condition cos0 + = — cos0 o {u + = u ). In this case, we have 

2u = w(T e /2). (35) 

Since w S i ne (T e /2) = 2/u e and Wtrian g (T e /2) = T e /4, we are led to the conclusion that the 
switching voltages are frequency independent and v ± me ' ) = ti£ nang ) = 1. 

As a way around this difficulty, we can apply a cutoff procedure where the equilibrium 
states are regarded as "saturated" states characterized by the cutoff parameter u s , tanh-u s m 
1. The result is that the relation between the angular and polarization parameters takes the 
modified form 

J tanhw, \u\ < u s , 

cos <p = < (3d J 

I sign(-u) tanh-u s = ± cos <p s , \u\ > u s 



and the parameter uq is now given by 



w(T e /2) T e < T s , 

^2u s = w(T s /2), T e > T s . V ; 



From Eq. (|3T|) it can be seen that the cutoff parameter u s determines the boundary frequency, 
f a = 1/T S , separating the regions of low and high frequencies. In the low frequency regime 
with f e = 1/T e < f s , the switching voltage parameters are given by 

v ( r e) = V(2 - u e /u a )u e /u a , <i trians) = y/JJTs, (38) 

where uo s = 2irf s . 

Note that, similar to the switching voltages, the switching time parameters, r + /T e and 
r_/T e , are frequency independent, r + /T e = 1 — r_/T e = 1/4, at f e > f s for all three wave 
forms (|4*])-(p)|). For the square- wave driving voltage of low frequency, these parameters assume 
linear dependence on the frequency: r + /T e = 1 — r_/T e = / e /(4/ s ). 



III. SWITCHING DYNAMICS 

From analysis presented in Sec. Ill Bl it can be inferred that properties of switching dy- 
namics are predominately determined by a certain class of periodic solutions (cycles) of the 
governing equation. Specifically, we have used the symmetry condition p + = —p-, where p + 
(p_) is the polarization parameter of the up (down) state, in combination with the cutoff 
procedure. 

Now we turn back to the model (j21j) and study how the anchoring energy influences the 
dynamics of switching. As in the preceding section, our attention will be focused on the 
frequency dependence of the switching voltages in the regime of complete switching. 

By contrast to Eq. (J23j) . the dynamical equation ({21} is generally not exactly solvable. 
Additional difficulties emerge in analysis of the solutions describing the switching process 
that need to be periodic in time and will be referred to as the cycles. As we shall be seeing 
the very existence of cycles does not follow from Eq. (J2TJ) immediately and the periodicity 
conditions may require somehow involved considerations. 

In this section we present an analytical approach to study the cycles based on approxi- 
mating the dynamical equation for the angular variable u defined in Eq. (|25j) . By applying 
this method, the case of square wave voltages can be treated analytically, whereas some 
relatively simple numerical analysis is needed for other wave forms. 

As a first step, we shall write the equation for the angular variable u 

du 

— = u — tanh(u) + [r(r) + w p ], p = cos0 = tanh(u) (39) 

deduced from Eq. (J2T|) by using Eq. (|23|) to make the change of variables: — > u. Then we 
apply the piecewise linear approximation for hyperbolic tangent 

tanh(u) - L(u) = h | U | < )' (40) 

I sign(Mj, \u\ > 1 

and obtain the approximate dynamical equation 

u = L{u) + [r(r) + w p ] (41) 

that will be the starting point for our subsequent considerations. Eq. (f4*TJ) retains the 
symmetry of Eq. (j3Uj) discussed at the end of Sec. Ill Al It can be readily checked that the 
transformation 

r — » r — T w /2, u — > —u, w p — > — w p (42) 



keeps both of these equations intact. 

It turned out that, as far as the cycles are concerned, the numerical results computed 
from Eq. ([39)1 are essentially the same as compared with the predictions of Eq. (|41|). In our 
calculations the relative error was found to be well below 0.1 per cent. 

Accuracy of the approximation is demonstrated in Fig. |3J It is seen that the difference be- 
tween the solid and dashed lines becomes noticeable only with a rise in numerical integration 
error as the interval of integration increases. 

As in Sec. Ill Bl suppose that the cell is initially in the down state with u(0) = —Uq 
and Uq > 1. After the first half-period time, the external voltage drives the cell into the 
up state with u{T w /2) = u + . This process might be called switching up and we consider 
it complete if u + > 1. By analogy, complete switching down takes place over the second 
half-period time when the cell goes from the up state with u(T w /2) = u + to the down state 
with u(T w ) = — m_ and U- > 1. 

Mathematically, this can be described in terms of the half-period mappings 

^ + :u = -u(0) ^u + = u(T w /2), 

^_ : u + = u{T w /2) |=> u_ = -u(T w ), (43) 

where and are determined by temporal evolution of the angular variable u over the 
first and second half-period, respectively. 

Owing to the symmetry ()42|) . \I> _ can be obtained from by changing the sing of 
the polar anchoring parameter w p . So, in the subsequent section, we concentrate on the 
half-period mapping 



A. Half-period mappings and cycles 

We restrict ourselves to the first half-period and study in detail how the angular parameter 
u evolves in time. For complete switching up, the solution can be written in the following 
general form 

(F (t), t<t , 
u{t) = \ F(t), r < r < n, (44) 
(F + (t), n<T<T w /2, 

where < r < T\ < T w /2 are defined by the conditions 

m(to) = -1, u(n) = 1. (45) 



It is not difficult to solve Eq. (j41)l and derive the explicit formulas for the functions that 
enter the solution (pl4*|) . The result is 

F (t) = -u + (w p -1)t + R (t), (46) 
F(t) = exp(r - 7o)K -l)-w p + exp(r)( J R(r ) - R(t)), (47) 
F + (t) = 1 + (w p + l)(r - 7i) + R (t) - Roin), (48) 



where 



Ro(t)= rir')dr f , R(r) = - / exp(-r') r(r')dr'. (49) 
Jo Jo 



From Eqs. (|44j) - ()48|) the half-period mapping from uq to u + can be derived in the 
parameterized form 



u -l=G (r ) = (w p -l)T + R (r ), (50a) 
u+ - 1 = G + ( n ) = (w p + l)(T tt /2 - n) + i2o(T tt /2) - i2o(n), (50b) 
£-( r o) = g±(r) = exp(-r)(w p ± 1) + R(r), (50c) 

where the parameters r and 7"i are defined by the two conditions (|45p. The first condition 
m(to) = -Fo( r o) — — 1 an d the relation u(T w /2) = F + (T w /2) = u + give u and u + as a function 
of r and Ti, respectively (see Eqs. ()5()aj) and (|50bjl for the corresponding expressions). 
The parameters To and t\ are coupled through the equation (|5()cj) which is an immediate 
consequence of the second condition from Eq. ()45j) : u[j\) = F(ti) = 1. 

Changing the sign of w p in Eqs. fl50aj) - (j50cj) gives the formulas for the half-period mapping 
Composition of the two mappings 



(51) 



relates the values of u at the end points of the period: «(0) = — u — ► u(T w ) = — u_. 

From the periodicity condition u = m_ it follows that, for cycles, % is the fixed point of 
So, we have the fixed point equations 

^(gi st) ) = <£ st) , <£ st) = «L st) - 1 = 4 st) - 1, (52a) 
^ + (gi st) ) =gJ t) = U J t) -l (52b) 

characterizing the up and down states of the cycle. The fixed point then can be found by 
solving Eq. (|52ajl . It provides the polarization parameter for the down state of the cycle: 
p_ = — tanhui* , whereas Eq. (J52bj) gives u+ and the polarization parameter of the up 
state: p+ = tanhw^ . 

Given the cycle Eqs. f)39 |l -()49 |l with uq = u_ yield u(t) for switching up during the first 
half-period, < r < T w /2. The reversal polarization current reaches a maximum at the 
instant of time r + where P(t+) = 0. 

For the second half-period time, T w /2 < r < T w , Eqs. (I39|) - (j49"|) with w = w+ and 
w p — > — w p give — u(t — T w /2). The current is now peaked at r = T w /2 + r_. From the 
condition p = and Eq. we can derive the equations for r + and r_ 

u(r±) [2u(r±) tanh M (r±) - l] « u(t ± ) [2u(r±)u(r±) - l] = r(r±), (53) 
u(t±) = u(r±) + [r(r±) ± w p ], (54) 

where < r-t < T w /2, u(r + ) = F(t + ) and u(r_) = F(t-)\ w -+- w ■ Then the switching 
voltages are given by 

r± = ±r(r±) = ±r e w ± , (55) 



where t>± = u(r±) is the switching voltage parameter. 
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FIG. 5: Critical period for square wave voltage as a function of (a) polar anchoring parameter, 
Wp = Wp/W and (b) voltage amplitude parameter, r e = PsV m /W. Thin dashed lines are computed 
from the approximate formula (|66|) . 



B. Square wave voltage 

Analytical results of the preceding section depend on the wave form through the functions 
defined in Eq. (JJDJ). The case of square wave voltages can be treated by replacing w p with 
w p + r e and setting R and R equal to zero. The result for the functions that enter the 
half-period mappings (see Eqs. (|50aJ) and (|50bJ) ) and the coupling equation ()50c|) is 

#± : G = (±w p + r e - 1)t , (56a) 

G + = (±w p + r e + l)(T w /2 - n), (56b) 

Ll fl a r e ±w p + l 

T\ = t + 7± = r + lnfj±, p ± = — -. (56c) 

r e ± Wp - 1 

By using Eqs. fl56aj) - ()56cj) along with the formulas (|5()ajl - (|50bjl the parameters To and n can 
be eliminated to yield the half-period mappings, ty + and as the linear functions given 
by 

q+ = V + (q ) = a+ - (3 + q = q[ +) (1 - q /q { +) ), (57a) 
g_ = vp_(g + ) = a _- (3_ q+ = q[-\l - q+fq^), (57b) 



where q^ = U; L — 1 and 

qf ] = a ± = (±w p + r e + l)(T w /2 - 7± ), q& = a ± /(3±. (58) 

Now that the expressions for the half-period mappings derived we discuss some of the 
most important consequences. The first point to address concerns the regime of complete 
switching. 

This regime occurs provided the parameters qo = Uq — I, q + = u + — I and g_ = u_ — 1 
are all positive. Since q = G (r ) with r > 0, from Eq. ()56a|) we obtain the inequality 

r e > 1 + \w p \ (59) 

as the complete switching condition for the driving voltage parameter, r e = PsV m /W. 

It is evident from Eqs. (|57a|) and (|57b|) that the parameters q + and g_ can be positive 
only if q[^ > 0. The latter and Eq. give the condition of complete switching for the 
period of driving voltage 



T w > T min = 2 max(7 + , 7_) = 2 In — — ^ (60) 



r e - 


w p \ 


+ 1 


r e ~ 


w p \ 


- 1 



restricting the driving frequency 1/T W from above. 

The conditions of complete switching (|59|) and (jBHj) define the complete switching thresh- 
olds for the amplitude and the period of the applied voltage. According to estimates that 
will discussed in Sec. IIV| the voltage parameter r e is typically much larger than unity and 
the asymptotic expression 

T min w 4(r e " 1 + KK 2 ) (61) 

can be used as a good approximation for the threshold period (ffiUj) at r e > 10. 

In general, the complete switching conditions are insufficient for cycles to exist. In other 
words, we need to impose more stringent constraints in order that the composition of map- 
pings (joTj) has fixed points. 

In our case, the mapping \1/ = o q — > q_ is linear and the parameter q varies in 
the range from zero to q^ . The fixed point exists only if the function q — ^(q) takes the 
values of opposite sign at the endpoints of the interval [0, q^]: q = and q = Since 

*(0) = [qt } ~ qt ] ]qV/qV, Q { o +) - ^(^ +) ) = ^ +) - q[~\ (62) 

the cycle condition \l/(0)[gg + ' ) — ^/(q^)] > can be conveniently written as the inequality 

P + P->0, P± = q^ -gi T) , (63) 

where the parameters q^ and q^ are defined in Eq. (fo'Hj) and linearly depend on T w . 

It can be shown that P + and P_ are both negative in the region where the cycles come 
into play. For this region of cycles, we have 

P± < at T w > T c , (64) 

where T c is the critical period given by 

C ~ (r e -l) 2 -^ +r ^ 1 1^1^ ln r 2_(i + hp | )2 - ( 65 ) 



Thus there are no cycles in the high frequency region above the critical frequency 1/T C . 

Similar to T min , the critical period T c can be conveniently approximated by the simplified 
asymptotic formula 

T c ^ 4(1 - I^D-V" 1 (66) 

in the region of high voltage parameters. Clearly, T c declines as the voltage amplitude 
increases, whereas the critical period becomes divergent when the magnitude of the polar 
anchoring parameter approaches unity. The curves depicted in Fig. El illustrate these effects 
and accuracy of the high voltage approximation. 

An important point is that, under the cycle condition (j64jl . the derivative of ^'(q) = 
[Qq Qo ] j is greater than unity. Therefore, the fixed points and the cycles are 
unstable. Note that stability of the cycles requires P + and P_ to be both positive. 

For the half-period mappings of the form (|57ajl and (|57b|) . it is straightforward to deduce 
the expressions for the parameters and q^ 

( St ) a± ~ f3±a T Qi^Qo^Pj ( ^ 

q± i-w- 1 ' 

that characterize the up and down states of the cycles. Solving Eq. (|53)l gives the switching 
time parameters r + and r + 

T ± = qf/(r e ±w p -l) + 5r ± , (68) 

= ln r e ±w p + ^(r e ± W y_ + 2 

2(r e ±w p -l) y 1 

that determine location of the reversal current peaks in time during the first and second 
half-period (switching up and down). The formulas 

u ( f « T w (l t w p ) [r e - (1 ± w p ) 2 r- x ] /4 + 2(1 ± 2 Wp )r e " 2 /3, (70) 
r± « T W {1 ± w p ) [1 + (1 T /4 T 4w p r e " 3 /3 (71) 

describe asymptotic behavior of w^ff ■* and r± in the region of high voltage parameters where 
r e > 10. 

As is expected from the symmetry 1)420. we have = u_ and r + = r_ in the limit of 
zero polar anchoring with w p = 0. This is also the case where the difference between the 
critical period T c and the threshold period T min disappears. 

Interestingly, at w p = 0, the leading term of the asymptotic expansion ([70)) [the first term 
in the square brackets on the right hand side of Eq. 1)70)1 ] can be derived from the symmetry 
condition (jHSj) discussed in Sec. Ill Bl More generally, Eqs. (ffH)) and (ffTjl indicate that, 
quite similar to the high frequency regime considered in Sec. Ill B) the parameters u± /T w 
and t±/T w are almost frequency independent having only a weak linear dependence on the 
frequency l/T w . 

The effects of the polar anchoring induced asymmetry depend on the polarity determined 
by the sign of w p . At the critical point with T w = T c and w p > [w p < 0], we have P+(T C ) = 

[P_(T C ) = 0] and the parameter q+ [q-^] vanishes. 

In the case where w p > and the down state is metastable in the absence of applied 
voltage, for the cycles, the polarization parameter of the down state appears to be 



higher than that of the up state p+ , p_ > p+ . An important consequence of this is that 
the switching time r + is longer than r_, r + > t_ and (r+ — tJ)/T w « u> p /2. So, at u> p > 0, 
the model predicts that, in cycles driven by the square wave voltage, switching down is 
faster than switching up. Note that the results for reversed sign of w p , can be obtained by 
interchanging the up and down states. 

C. Sine- wave voltage: numerical analysis 

Now we extend the analysis of the previous section to the cases of sine- wave and triangular 
voltages. More specifically, we shall present the results for the case of the sine-wave voltage. 
The triangular driving voltage can be treated in just the same way but this involves rather 
cumbersome expressions. So, our subsequent considerations are equally applicable to both 
of the waveforms. 

As a first step, we perform integrals in Eq. for the waveform function (J3J) and deduce 
the expressions for the functions 



R{Q = r e exp(-£/a; w )[a; tt cosf + sinf]/(l + (72) 

#o(0 = r e[l - cos£\/uj w , £ = uj w t, (73) 

where u w = 2n/T w = utw is the renormalized frequency (frequency parameter). Substi- 
tuting Eqs. (|72*|) and ([75]) into Eqs. ()50a|) - ()50cj) gives the half-period mapping in the 
following parameterized form 

^+ : q = u w {u - 1) = (w p - 1)6 + r e [l - cos£o] = Gf^fo), (74a) 

q + = cv w (u+ - 1) = (w p + 1)(tt - 6) + r e [l + cos 6] = G^fo), (74b) 

y-(Co) = </ + (e0, fc(0 = exp(-e/^)K ± 1) + i2(0, (74c) 



where the parameters ^ = u w t and ^ = uj w ti are ranged between zero and tt. [Note that 
the parameters g« differ from those defined in Sec. IHI Bl bv the factor u w .} 

As opposed to the case of square wave voltage, the parameters £o an d Ci cannot be 
eliminated. We begin with locating the values of these parameters that represent the regime 
of complete switching where qo and q + are non-negative, go,+ ^ 0. 

Our first remark concerns behavior of the functions, Go,+ and g±, that enter the mapping 
^> + . The time derivatives of these functions 

-l±K + r(r)), (75) 

exp(-r)(tw p ±l + r(r)) (76) 

show that over the first half-period time G + and g + are monotonically decreasing functions 
of r. 

By contrast, the functions Go and — g_ exhibit nonmonotonic behavior. As r increases, 
they first decay reaching the minimum located at r min and then grow up to the maximum 
at t = r max decreasing at r > r max . The points r min max are determined by the equation 
G (r) = [or, equivalently, v(r) = (1 - w p )/r e ] and G (r min ) < 0. 



<9G , + (r) 
Or 

dgjj) 

dr 



For the sine-wave voltage, £ min = arcsin[(l — w p )/r e ] and £ max = 7r — £ min . The function 
Gg™ e ' monotonically increases from G sme ^ (£ m i n ) < to G sme ' ) (£ max ) and it passes through 
zero, Go Sint ^(£o°^) = 0> before the waveform function reaches a maximum onlyifG[, sine) (7r/2) > 
0. The latter provides the complete switching condition for the driving voltage parameter 

r e > (l + K|)7r/2. (77) 

Interestingly, replacing r e with nr e /2 in Eq. (J77|) recovers the result for the square wave 
voltage given by Eq. (JHSJ). 

The largest value of £o, £o = £,o + \ can be found by solving the coupling equation (|74cj) 
at £i = 7i. Owing to the inequality l;^ > £q°\ the complete switching condition for the 
frequency is given by 

gW) > g+{*). (78) 

It can be checked that Eq. (J5UJ) with r e replaced by 2r e / n provides a good estimate for T min 
defined by the condition (|7Sjl. 

(0) ^ t ^ t(+) 



At this stage we have found that the endpoints of the interval for £o : £o < £o < 6 



o > 



can be evaluated as solutions of the equations 

9-ffl)=9 + W> Gt c \^) = 0. (79) 

According to the coupling equation (|74c|) . when £ increases from ^ to £q + \ the parameter 
£i changes from £^ to ix. The equation for is 

9 + (£ +) ) = 9-{$ ) )- (80) 

For the parameters go and q + , it can be inferred from the above results that q + mono- 
tonically decays from q[ + ^ to zero with q growing from zero to q^ . The values of q^ and 
q[ + ^ can be computed from the relations 

g (+) = q w = c^)^W )t (si) 

where £q + ' ) and ^[ + ^ are given in Eq. (J7HJ) and Eq. (jSOjl. respectively. 

So, we arrive at the conclusion that the parameterized mapping (|74|) behaves just like the 
half-period mapping for the square wave voltage (|57a|) . In addition, numerical calculations 
show that using the formulas (|57|) as a linear approximation for the parameterized mappings 
ty + and \tf_ does not introduce any noticeable errors. [The relative error was typically below 
10 V 

It immediately follows that the cycle condition (jfiHJ) and Eq. (Jfi7j) remain applicable for 
the sine-wave voltages. However, the switching time parameters £ + and £_ can only be 
found numerically by solving Eq. (JHHj) where the relations (fT2~J) and ()74aj) are used to define 
the function (|4"7|) and to evaluate the values of £o for q± l \ respectively. 

For the most part, qualitative predictions of the model for the sine- wave voltage remain 
unaltered as compared to the case of the square wave voltage. The cycles are found to be 
unstable and the condition (|64|) determines the cycle region. Despite the critical period 
is longer than T c given by Eq. (|65|h its dependence on the voltage and polar anchoring 
parameters is qualitatively the same. 



Similar to Eq. (|70J). the T- dependence of the parameters u± turned out to be approx- 
imately linear, whereas the switching time parameters £± show a weak dependence on the 
frequency. The effects of polarity are analogous to those discussed at the end of Sec. IIII Bl 

Thus, for all the wave forms under consideration, the fixed points of the mapping ()51|) are 
repelling and the resulting instability of cycles is the characteristic feature of the model (J39|) . 
Strictly speaking, it means that initially small deviations from a cycle will grow in time. 

In Appendix |X] our analysis is extended to the generalized model with the double- well 
potential which is a continuous deformation of Eq. Equivalently, it can be regarded 

as a modification of the Rapini-Papoular potential. The key results can be summarized as 
follows 

(a) There are two limiting frequencies, f\ and f 2 , so that the cycles are unstable at f 2 < 
f w < fi (see Eq. flA~18jl and Eq. (jMj ). 

(b) There are no unstable cycles until the separation between the minima of the potential 
is sufficiently large (see the condition ()A16j) ). 

(c) In addition to unstable cycle, there is the branch of stable cycles in the low frequency 
region below the critical frequency f c , f w < f c (see Eqs. (IA20j) - (jA21|) h 

(d) At w p 7^ 0, the low frequency region is separated from the frequency interval of unstable 
cycles by the gap of the width f 2 — f c . 

These results suggest that the Rapini-Papoular anchoring potential breaks down at large 
angular deviations of the director rendering the cycles unstable. For similar reasons, mod- 
ifications of the potential were considered in very recent studies dealing with the dynamics 
of a pitch jump in cholesteric liquid crystal cells |32l |33|. 
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Tilt angle 


Viscosity 


Mixture name 




Phase sequence 




(nC/cm 2 ) 


9 (deg) 


7 (Pa s) 


FLC-497 




Sm-C*^> Sm-A*^* 


Iso 


95 


27 


0.11 


FLC-510 


Cr ^ 


Sm-C*^* Sm-A*^* 


Iso 


98 


31 


0.18 



TABLE I: Parameters of the FLC mixtures measured in the Sm-C* phase at the temperature 23°C. 



IV. EXPERIMENT 

In this section we present the experimental results on frequency dependence of the thresh- 
old voltages measured in asymmetric SSFLC cells. It is found that the asymmetry induced 
shift of the hysteresis loop increases with the driving frequency. We apply the theoretical 
results of Sec. IIHI to model the process of switching within the cells and to estimate the 
polar anchoring parameter from the experimental data. 



A. Sample preparation 



In our experiments we used asymmetric FLC cells where the FLC layer is sandwiched 
between two dissimilar substrates. One of the substrates was covered with a photo-aligning 
substance - the azobenzene sulfric dye SD-1, whereas the other one was simply washed 
in N,N-dimethylformamide (DMF) and covered with calibrated spacers. By contrast, the 
substrates with photoaligned films identical in anchoring properties were assembled to form 
symmetric FLC cells. 

Following the procedure described in Ref. jsij, SD-1 was synthesized from corresponding 
benzidinedisulfonic acid using azo coupling. The solution was spin-coated onto glass sub- 
strates with indium-tin-oxide (ITO) electrodes at 800 rpm for 5 seconds and, subsequently, 
at 3000 rpm for 30 seconds. The solvent was evaporated on a hot plate at 140°C for 10 
minutes. 

The surface of the coated film was illuminated with linearly polarized UV light using a 
super-high-pressure Hg lamp through an interference filter at the wavelength 365 nm and 
a polarizing filter. The intensity of light irradiated normally on the film surface during 30 
minutes was 6 mW/cm 2 . 

We used two different pitch-compensated liquid crystal mixtures: FLC-497 and FLC-510 
(from P. N. Lebedev Physical Institute of Russian Academy of Sciences) as materials for 
the FLC layer. The mixtures were injected into the cells in an isotropic phase by capillary 
action at T=85°C and T=100°C for FLC-497 and FLC-510, respectively. In these mixtures, 
the FLC helix is unwound due to compensation by two chiral dopants with opposite sense of 
chirality (opposite sign of handedness) but the same sign of the spontaneous polarization . 
The parameters of the mixtures are listed in Table HI 



B. Experimental results and modelling 

Measurements of the hysteresis loops were performed using the standard electro-optical 
set-up composed of a He-Ne laser, a Hewlett Packard Infinum digital oscilloscope and a 
generator of triangular pulses. We also used a rotating table for adjusting the angular 
position of the FLC cell placed between crossed polarizers. The voltage amplitude of the 
generator, V m , used in all experiments was 10 V and the driving voltage frequency can be 
varied in the range from 10~ 4 to 10 3 Hz. 

Figure |3] schematically represents the form of a typical hysteresis loop measured in our 
experiments. The threshold voltages, V±, and the corresponding values of the switching 
voltage parameters, v± = ±V±/V m , then can be extracted from the experimental data for 
the electro-optical hysteresis loops. For asymmetric cells filled with FLC-510 and FLC-497, 
the results shown in Fig. El clearly indicate the difference between the voltage parameters 
for switching up and down, v + and v_. Therefore, the hysteresis loops are shifted. By 
contrast, for symmetric cells, the shift, v + — t>_, is found to be vanishing within the limits 
of experimental error. 

The curves presented in Fig. |3 were measured in the cells where the thickness of the 
FLC layer, d, and of the aligning film, d^ was 5 /xm and 12 nm, respectively. In our 
previous papers ^M^, E3, E^t , we studied anchoring properties of the azo-dye films and our 
experimental technique can be applied to measure the strength of anchoring W in FLC cells. 
It was estimated to be about 2 x 10 -4 J m~ 2 and 1.2 x 10~ 4 J m -2 in the cells filled with 
FLC-510 and FLC-497, respectively. 



0.8 



0.6 



<U 
SJO 

eg 

*?0.4 

s 

in 



0.2 



~ i — i 1 1 ii 1 — i 1 1 in 1 — i Mill 1 — i Mm 1 — | p 1 — i — i — i r 

FLC-510 

A 

A A V+ 

T T V_ i 



(a) 



J I I I I L 



uL 



.01 0.1 1 10 100 

Frequency (Hz) 



I I I I I I I I I I 

FLC-497 


i i i i i | 






A 


A A 


v+ 


A _ 


T T 


v_ 


A 

a' - 


- 










a 7 * t 


(D) 




If 




A 

i 






a' 




- 




W — 

▼ 






r 




















i i 




1 



0.9 



0.8 



0.7 > 

ii 

0.6 M 
"o 

0.5 M 
S 

u 

0.4 

CO 



0.3 



10 



100 



0.2 



Frequency (Hz) 



FIG. 6: Experimental curves for switching voltage parameters, v± = ±V±/V m , measured as a 
function of frequency in the cells filled with (a) FLC-510 and (b) FLC-497. 



The characteristic times tg and tw given in Eq. (|2U|) can now be estimated as follows: 
t E ~ 0.9 x lCr 4 s and t w « 4.5 x 10~ 3 s for the FLC-510 cell; t E » 0.58 x 10" 4 s and 
tw ~ 4.58 x 10 -3 s for the FLC-497 cell. The corresponding values of the driving voltage 
parameter are: r e « 50 and r e w 79. 

The above estimates define the parameters that enter the dynamical equation (J2TJ) of the 
model discussed in Sec. Ill Al In this model the director field (0) is assumed to be uniform 
across the cell. Such assumption can be a reasonable approximation in the high field regime 
of switching. This is the case in which the characteristic length = (K/PsEq) 1 ^ 2 , where 
K is the effective elastic constant, is shorter than the cell thickness, £e < d. From the 
estimate m 0.07 /im <rf = 5.0 fim obtained using a typical value of the elastic constant, 
K Pd 10~ n J m^ 1 , we may conclude that our measurements were performed in the high field 
regime and the uniform model can be used to interpret the experimental results. 

We can also estimate the anchoring extrapolation length £w = K/W ~ 0.05 — 0.08 fim, 
so that the characteristic lengths £e and £w are of the same order. Therefore, the anchor- 
ing conditions at the boundary surfaces described by the anchoring energy could play an 
important part in the switching dynamics. 

In Sec. Ill A| we found the expression for the electric field inside the FLC layer (|T5|) 
which differ from the external electric field due to the presence of the insulating aligning 
film. This effect crucially depends on the value of the dimensionless parameter r] defined in 
Eq. (fT3j) that, in our case, can be estimated by assuming that d<i = 0, d\ = dA = 12 nm and 
E\ = sa ~ 7.5 is the dielectric constant of the SD-1 layer. For both FLC-510 and FLC-497 
cells, we obtain small values of rj that are below 5 x 10~ 4 . As a result, even a rough estimate 




Frequency tw/T Frequency tw/T 

FIG. 7: Frequency dependence of parameters (a) — 1)/T W and (b) (u+ — T)/T w computed from 
the experimental data for the FLC-510 cell using the model (|41|) with r e = 50 and w p = —0.65. (The 
parameter u± determines the polarization of the up (down) state P± = ±Ps tanhu±). Piecewise 
linear approximation of the results (solid lines) defines the ranges of high and low frequencies with 
the boundary frequency about 0.07 /tw 16 Hz. 

for the depolarizing voltage Vd ep ~ -PsgU/^a^o) ~ 0.16 V derived from Eq. (|THj) yields the 
voltages an order of magnitude lower than typical values of the voltage shift measured in 
our experiments. In addition, a more accurate analysis of the contributions to the effective 
electrostatic free energy (JTfiJl clearly shows that the voltage drop across the aligning layer 
cannot be responsible for asymmetry of the switching voltage thresholds, V + ^ V— 

The method developed in Sec. IHII can now be applied to model switching dynamics of 
the FLC cells using the experimental results. To this end we calculated the parameters 
u± characterizing the polarization of the up and down states from the data on frequency 
dependence of the switching voltage parameters v±. This procedure involves solving Eq. (JHHJ) 
and using the half-period mapping (|5()jl for the triangular waveform function (jSJ. The 
mapping \1> + ( 1 i r _) is used to compute the parameters u± for switching up (down) from the 
data on the switching voltage v + (i>_). 

In general, the computed values of u± for switching up deviate from the corresponding 
results for switching down. We used the polar anchoring parameter fitting parameter 

to minimize the difference and to obtain the results nearest to the cycle. For the FLC-510 
cell, the value of w p is found to be about —0.65 and we present the results for the parameters 
(u± — l)f w in Fig. Referring to Fig. there are two frequency regions where, similar to 
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FIG. 8: Switching voltage parameter computed as a function of frequency (solid lines) using linearly 
approximated frequency dependence of (u± — l)f w for (a) the FLC-510 cell and (b) the FLC-497 
cell. 

cycles, the frequency dependence of (u± — l)f w is nearly linear. Using linear approximation 
for the parameters (u± — l)f w we estimated the boundary frequency separating the low and 
high frequency regions at about 16 Hz and calculated the curves for the switching voltage 
parameters shown in Fig. as solid lines. 

The switching voltage vs frequency curves computed along the same lines for the FLC-497 
cell are given in Fig. Bp. In this case, the polar anchoring parameter is found to be about 
—0.78. From Fig. it can be seen that the linear approximation in the high frequency 
region reproduces the experimental data for the voltage coercitivity remarkably well. The 
scatter in the data and its scale are more evident in the results for the voltage shift shown 
in Fig. 

Despite good agreement between the theoretical and the experimental results, it is es- 
sential that the linearly fitted data, strictly speaking, do not represent the cycle and the 
modelling was only performed over the first period. Since the cycles are unstable, evo- 
lution in time will have the deterioration effect on the agreement between modelling and 
measurements. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied how the dynamical characteristics of switching in SSFLC 
cells are affected by asymmetry of dissimilar substrates in anchoring properties. The asym- 
metry effects are found to be dominated by the polar contribution to the anchoring energy 




Frequency t w /T Frequency t w /T 

FIG. 9: (a) Voltage coercitivity and (b) shift of the hysteresis loop as a function of frequency for 
the FLC-497 cell. Solid lines represent the results of modelling at r e = 80 and w p = —0.78 using 
the approximate relations: — l)f w ~ 0.15 + 4.7 f w and (u + — T)f w ~ 15.6 — 4.24f w . 

and manifest themselves in the frequency dependent shift of the hysteresis loop. 

Our theoretical considerations were primarily concerned with the predictions of the uni- 
form theory on the frequency dependence of the switching voltage thresholds (the switching 
voltage parameters). It was assumed that the steady state regime of switching is determined 
by the cycles which are time-periodic solutions to the dynamical equation. So, we have de- 
veloped the method to analyze the properties of cycles depending on the driving frequency. 
In this method the dynamical equation for the angular variable is approximated to yield 
the half-period mappings in the parameterized form. The cycles then studied in terms of 
the fixed points of the composition of two half-period mappings. 

By using this method we performed analysis for the cases of square wave, sine-wave and 
triangular voltages. It was found that the cycles are unstable and can only be formed when 
the driving frequency is lower than its critical value, f < f c — 1/T C . The critical frequency 
f c declines with the polar anchoring parameter w p suppressing the cycles at \w p \ = 1. 

The polar anchoring breaks the mirror symmetry relating the processes of switching up 
and down leading to the difference in the magnitude of the corresponding switching time 
and switching voltage parameters. So, we have r + ^ r_ and v + ^ v_ at non-vanishing 
polar anchoring parameter, w p ^ 0. Therefore, the shift of the hysteresis loop, v + — i>_, 
results from the symmetry breaking effect caused by the polar contribution to the anchoring 
potential (JHJ) characterized by the polar anchoring parameter ©. 

Our calculations indicate that, for cycles, dependence of the parameters u±/T and t±/T 



on the frequency is nearly linear. Modelling of the switching dynamics using the experimental 
data revealed similar behavior provided the frequency is not too low. Since the conductivity 
of the FLC layer has been neglected in the electrostatic model discussed in Sec. Ill A[ we 
cannot expect the model to give accurate results at low frequencies. 

The expression for the effective electrostatic potential (|TH|) shows that there are no sym- 
metry breaking contributions due to the voltage drop across the insulating aligning films 
and the dielectric anisotropy of the FLC layer. So, they cannot be responsible for the asym- 
metry effects under consideration. Thus we may conclude that the anchoring energy is the 
determining factor in this problem. Of interest is the fact that similar conclusion underlines 
the method proposed in [39( to control the process of switching. 

During the switching process the azimuthal angle varies between zero and 7T. This involves 
large deviations of the director from the easy axis, so that the switching dynamics can be 
sensitive to the shape of the anchoring potential. It can be assumed that the instability 
of cycles is a characteristic of the model with the potential taken in the standard Rapini- 
Papoular form ji^ . In the Appendix we have studied the generalized model with the double- 
well potential and found that the branch of stable cycles emerges even at small deformations 
of the potential. Note that the results derived in the Appendix not only demonstrate the 
effectiveness of the analytical approach but also can be placed in a more general physical 
context. 

The uniform theory cannot be directly applied to more complicated cases with spatially 
inhomogeneous orientational structures involved. So, discussion of such problems as the 
soliton (kink) mechanism of switching 41, 4^, H| and dynamics in the chevron geometry [Zil 



1451 ] is beyond the scope of this paper. Nevertheless, it is pertinent to note that, under 
certain circumstances, the switching dynamics of inhomogeneous structures can be effectively 
reduced to a uniform model for relevant spatially independent variables using a trial solution 
in the first approximation. 
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APPENDIX A: DOUBLE- WELL POTENTIAL: CYCLES IN PIECEWISE 

QUADRATIC MODEL 

In Sec. IHII we have studied the switching dynamics by analyzing the dynamical model (}4T| 
for the angular variable u. The assumption underlying the study is that the dynamics 
is governed by the cycles which are periodic solutions to the dynamical equation in the 
regime of complete switching. Technically, we have used the analytical approach described 
in Sec. IIII Al in which the cycles are related to the fixed points of the composition of the 
half-period mappings (JoTj) . 

One of the most important results presented in Sec. IIII Bl and Sec. IIII CI is the instability 
of the fixed points, so that there are no stable cycles in the model (j41j) . In this appendix, in 
order to clarify the mechanism of such instability, we consider the generalized model with 
the governing equation 

du dU{u) , , TT , 

q^ = u = ^T +r ( r )' U(u) = U (u) - w p u, (Al) 



where the anchoring energy potential U is taken to be a piecewise quadratic function of the 
following form 

(e„u 2 /2 + (1 + e_)(u+ 1/2), u<-l, 
U (u) = < -u 2 /2, -1 < u < 1, (A2) 

[e+M 2 /2- (l+ e+ )(u- 1/2), u > 1. 

Clearly, our original model (jH| can be regarded as the limiting case of the modified 
model (jAlj) in which the parameters e± are equal to zero, e± = 0. 
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FIG. 10: Double- well effective potential as a function of u. Two cases are shown: (a) e± = 1 at 
various values of the polar anchoring parameter w p ; and (b) e+ 7^ e_ at u> p = 0. 

When e± > 0, [/ is a double- well potential that have two local minima representing 
the equilibrium up and down states with u = u^ n and u = u^ n , respectively. These are 
separated by the maximum located at u = « max and are given by 

= ±(1 + e± ± w p )/e±, 2U^l = 1 + e ± - (1 + e± ± w P ) 2 A±, (A3) 

"max = "Wp, ^max = U>J/2. (A4) 

From Eq. ()A3|) the minima are energetically equivalent and the states characterized by u^ n 
and u^ n = — are stable at e + = e_ and w p = 0. This is no longer the case when w p 7^ 0. 
As is illustrated in Fig. ITUk . similar to the anchoring potential (jBJ) [see Fig. |2J, the down 
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FIG. 11: (a) The functions P± versus the period T w computed at w p = 0.2, r e = 50 and e = 
0.5 < e c ~ 0.8. (b) The critical value of e as a function of the voltage parameter at various 
values of the polar anchoring parameter w p . The dashed curves are computed from the asymptotic 
formula (IA15I). 



state is metastable at w p > and the energy barrier disappears when the polar anchoring 
parameter approaches unity. 

Referring to Fig.llOb. the asymmetry effects can also be caused by the difference between 
e+ and e_. But, for brevity, we shall restrict ourselves to the case specified by the conditions 



e± = e > 0, \w p \ < 1. 



(A5) 



In the limit of vanishing e, e — > 0, the parameters become infinitely large and 

the minima shift at infinity in opposite directions. We may anticipate that this divergence 
renders the cycles unstable and our study of the model with the double-well potential (|A2j) 
will confirm this hypothesis. 

We perform our analysis following the general procedure described in Sec. IIII Al For the 
dynamical equation (jAljl . the structure of the solutions is defined by Eq. (jUJ), where the 
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FIG. 12: (a) Region of unstable cycles, fa < f w < f%, in the w p -f w plane. At r e = 5.0 and 
r e = 10.0, the curves for the limiting frequencies, fa and fa, are shown as solid and dashed lines, 
respectively, (b) The region, fa < fw < fa, in the r e -f w plane. Solid and dashed lines represent 
the curves for the limiting frequencies computed at w p = 0.2 and w p = 0.4, respectively. 



functions Fq and F + need to be changed as follows 

F (t) = — u exp[— er] + (w p — 1 — e)(l — exp[— er])/e + exp[— er] R e (r), (A6a) 
F+(t) = exp[-e(r - n)] + (w p + 1 + e) (l - exp[-e(r - n)])/e+ 

exp[-er] (^(r) - Rein)), (A6b) 

7?,(t-)=/ exp(er')r(r')dr. (A6c) 



The corresponding modification of the functions Go and G + that enter the half-period map- 
ping (|5U|) is given by 



G (t ) =(w p - l)(exp[er ] - l)/e + R e (r ), 
G+(n) =exp[-eT w /2]{(w p + 1) (exp[eT w /2] - exp[en])/e+ 
R e {T w /2)-R e {n) 



(A7a) 



(A7b) 



Note that the coupling equation (|50cj) as well as the function R denned in Eq. (149)1 remain 
untouched. 
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FIG. 13: Critical frequency of stable cycles, f c = 1/t] s) [see Eq. jS2DJ], 

as a function of (a) the 

polar anchoring parameter and (b) the parameter e. The asymptotic formula (|A21|) is used for 
computing the curves shown as thin dashed lines. 



1. Square wave voltage 
For the square wave voltage, Eqs. (|A7a|) and (|A7b|) take the simplified form 

G (r ) =(w p + r e - l)(exp[er ] - l)/e, (A8a) 
G+(n) ={w p + r e + 1)(1 - exp[e(r 1 - T w /2)])/e. (A8b) 

These expressions can now be combined with the coupling equation ()56c|) to yield the linear 
half-period mapping of the form (|57jl. But, in place of Eq. (J55j) . we have 

gS ±} = (±w p + r e + 1) (1 - exp[e( 7± - T w /2)])/e, (A9a) 
^ = (±u> p + r e - l)(exp[e(T w /2 - 7± )] - l)/e, (A9b) 

where j± is given in Eq. (|5ficjl . 

Following the same line of reasoning as in Sec. IIII B| we derive the conditions of complete 
switching for the voltage parameter and the frequency (see Eq. (JHHj) and Eq. (jEOJ), respec- 
tively). We also find that the cycles may only develop provided the cycle condition (jfiHJ) is 
met. 
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FIG. 14: Frequency dependence of parameters u± /T w for (a) stable and (b) unstable cycles with 



w p = 0.2 and r e = 50 at various values of e. 



This condition is expressed in terms of the parameters P±. Substituting the relations (|A9 
into Eq. ()63|) yields these parameters 



ezP± = a±z 2 — 2r e z + b±, z = exp(— eT w /2), 

a± = (Tw p + r e + l)[(3 T }\ 6± = (±w p + r e - 1)\J3±] 



(A10) 
(All) 



where f3± is defined in Eq. IjSficJ) . It remains to recall that the cycles may build up only if 
P+P_ > 0. 

From Eq. ()A10j) it is clear that the sign of the parameters P± is determined by the 
quadratic functions of the period dependent variable z that lies in the interval 



< z < z- 



niax; *max 



exp(-eT min /2), 



(A12) 



where T min is given by Eq. (foTIj) . 

The discriminants of the quadratic polynomials (|A10|) 



D ± = r 2 e -(r 2 e -(l T w p ) 2 ) 



e 



rj — (1 ± w p ) 2 



are both positive when the parameter e is below its critical value given by 
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FIG. 15: Frequency dependence of switching time parameters t±/T w for (a) stable and (b) unstable 
cycles with w p = 0.2 and r e = 50 at various values of e. 

In the region of high voltage parameters, the asymptotic formula for e c is 

e c « \w p \-\l - \w p \) 2 (2 - (1 + K|) 2 r; 2 )/8. (A15) 

Now, assuming, for definiteness, that w p > 0, let us see what happen if the condition 

e < e c (A16) 

is fulfilled. It can be shown that P+(T min ) = P + (z max ) > and P_(T min ) = P-(z max ) < 0. 
Referring to Fig. fTTk. P + changes its sign as T w passes through the point and is negative 
on the interval (T± u \ T^). The endpoints of this interval can be found as the roots of the 
polynomial ezP + (z) and determine the region of unstable cycles 

P ± < at T[ u) <T w < T 2 (n) , (A17) 

where 

T 4 (n) = -2c" 1 In*, z h2 = (r e ± y/D^)/o+ (A18) 

and P+(*) = 0. 

Outside the interval (|A17[) . P + is positive and the root of ezP_(z), P-(z c ) = 0, gives the 
point T w = Tc^ where the sign of P_ is reversed. So, in addition to the unstable cycles, we 



have the low frequency region of stable cycles 



P± > at T w > T c (s) , (A19) 

where 

T (s) = _ 2e - 1 lnz c , z c = (r e - ^Dl)/a^. (A20) 

The asymptotic formulas for the critical period of stable cycles T^f 1 and the limiting periods 
of unstable cycles T{1 are 

T« « 2e~ 1 r~ 1 (l + \w p \ + 2e + ^/(l + \w p \) 2 + 4e\w p \) , (A21) 
T$ « 2e- 1 r e - 1 (l - | Wp | + 2e T ^/(l - Kl) 2 - 4e|w p |) . (A22) 

The corresponding frequencies, f c = 1/Tc and f t = might be called the critical 

frequency of stable cycles and the limiting frequencies of unstable cycles, respectively. These 
frequencies define two regions: (a) the low frequency region of stable cycles with f w = 
1/T W < f c and (b) the high frequency region of unstable cycles with f 2 < f w < ft. 

The limiting frequencies, fi and f 2 , merge together at e = e c and the branch of unstable 
cycles disappear as the parameter e further increases. Fig.lllb shows that the critical value of 
e, e c , increases with the voltage parameter r e . By contrast, as it can be seen from Eq. (jA15|) . 
the polar anchoring reduces the value of e c thus suppressing the unstable cycles. 

So, there are no unstable cycles if the polar anchoring parameter is sufficiently large. 
The inhibition of unstable cycles by the polar anchoring is clearly evident from the curves 
presented in Fig. H21 

Similar to the limiting frequency f\, the critical frequency f c declines as the parameter 
w p increases (see Fig. IT3k). But the frequency f c does not vanish and the stable cycles are 
not suppressed at \w p \ = 1. 

In the limit of zero polar anchoring with w p = 0, the limiting frequencies, fi and f 2 , 
are equal to the threshold frequency / min = 1/T min , f\ = f min , and the critical frequency 
f c , correspondingly. The result is that the stable and unstable cycles are separated by the 
branching point f 2 = f c - 

Splitting of the branching point occurs as the effect of the polar anchoring induced asym- 
metry when w p ^ and f 2 > f c . It means that there is the gap separating the frequency 
intervals for stable and unstable cycles. The width of this gap, which is the frequency 
difference f 2 — f c , grows with the polar anchoring parameter. Interestingly, similar effects 
take place for wave numbers of helical structures in asymmetric cholesteric liquid crystal 
cells |46[. 

From the above discussion it follows that, when the parameter e increases, the limiting 
frequencies move towards each other until they coalesce at the critical point e = e c . By 
contrast, the critical frequency f c is an increasing function of e (see Fig. 113b ) and the 
frequency range for the stable cycles expands with the parameter e. 

In the opposite case where e — > 0, the frequencies f c and f 2 are both vanishing along with 
the branch of stable cycles. In accord with the results of Sec. IIIIB| the limiting frequency 
fi approaches the critical frequency l/T c (see Eq. fjfiH]l ) for the model (JHJ). 

The parameters q± = — 1 characterizing the up and down states of the cycles can be 
computed by substituting Eq. (jA9a|) and Eq. (jA9b|) into the relation (|67p. For the switching 



time parameters, the first term on the right hand side of Eq. (|68|) needs to be changed by 
using the function (|A8a|) . The result is 

r± = e" 1 ln[l + e q ( f/(r e ±w p - 1)] + Sr ± , (A23) 

where <5r± is given by Eq. (|69jl. It can be shown that the asymptotic expressions for r± 

er ± « ln[2/(l + z)] + ^((1 + e)(l - z)/(l + z) T w p ), (A24) 

where z is defined in Eq. (IA10)) . are applicable only for stable cycles at sufficiently low 
frequencies. 

Dependence of the parameter u± /T w on the frequency is shown in Fig.lT4"k and Fig. IT^b 
for stable and unstable cycles, respectively. Similar results for the switching time parameters 
t±/T w are given in Fig. O 

It is seen that the curves for the stable and the unstable branches reveal qualitatively 
different behavior. In particular, by contrast to the unstable cycles, for the stable cycles, the 
polarization of the up state is higher than that of the down state and switching up is faster 
than switching down. This effect can be loosely described as the cycle is shifted towards 
either stable or metastable state depending on the cycle stability. 

Our concluding remark concerns the sine-wave and triangular voltages. These cases can 
be studied in just the same way as for the model (j41J) using the procedure described in 
Sec. IIII CI A comparison between the time derivatives of the modified functions 



^Mll = ( Wp -l + r (r)) exp[6r], (A25) 

= -(w p + 1 + r ( T )) exp[e(r - T w /2)\ (A26) 

and Eq. (fTH|) shows that the change of the analytical expressions does not affect behavior of 
the functions and the general results of Sec. IIII CI are equally applicable to the model with 
the double- well potential ()A2j) . After performing numerical analysis for these waveforms, 
we have found that the predictions of the model discussed in this appendix have not been 
changed qualitatively. 
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